Interband Effects of Magnetic Field on Hall Effects for Dirac Electrons in Bismuth 
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Hall effects are investigated for three-dimensional Dirac electrons as a model of bismuth alloys. 
It is found that there is unconventional contributions to the Hall conductivity a^y generated by 
the interband effects of a magnetic field. This phenomena is remarkable near the band-edge. The 
Hall coefficient exhibits two unexpected properties; a sharp peak at around the band-edges, and a 
drastic sign change in the band gap region. Implications of the present results to bismuth alloys are 
discussed. 
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Electrons in a periodic potential can be described as 
Bloch band theory, which has been proved very useful for 
many purposes. In this theory, it is natural to expect that 
the effects of an external magnetic field H are described 
by the effective Hamiltonian ,3^cS = En{k + eA/c) ex- 
cept for the Zeeman effect, where the ?i-th Bloch band 
is characterized by an energy function En{k) for H ~Q. 
This simple procedure has been proved to be very use- 
ful for interpreting many phenomena. Nevertheless, this 
Bloch band picture is insufficient in principle[l|, |3]; elec- 
trons under a magnetic field are not confined in a single 
Bloch band, but undergo complex interband oscillations. 
These effects will be referred to as interband effects of 
a magnetic field, which are not contained in the Bloch 
band picture. Although the interband effect is a fun- 
damental problem even in the non-interacting electron 
systems, it has been studied only in limited cases [3[, i.e., 
the orbital magnetism, where the interband effect plays 
a crucial role^j. In this letter, we study the interband 
effects on the Hall effects having Bi and its alloys in mind 
as a typical experimental situation. 

Bismuth and its alloys have been famous for their 
large diamagnetism for a long time. Especially they take 
largest values when the chemical potential is located in 
the region where the density of states (DOS) vanishes[5|. 
This can never be understood in the Bloch band picture, 
i.e., the Landau-Peierls formula|6|, which indicates that 
the orbital susceptibility x is proportional to the DOS. 
This mystery has been solved by taking account of the 
interband effect [J|. In particular, Bi is a semimetal and 
its band structure has a special form, i.e., a narrow-gap 
Dirac electron systems (Fig. [T]). This is the reason why 
the interband effects on orbital magnetism are prominent 
inBi. 

The Hall conductivity, Uxy, is more complex than x- 
The orbital susceptibility originates from the permanent 
orbital current, which is non-dissipative. On the other 
hand, Uxy in weak magnetic field is essentially dissipative. 
In the following, we study cixy in a three-dimensional 
Dirac system to show that Cxy is actually related to the 



diamagnctic current. Our results give some clues to a 
long-standing problem about how electrons contribute 
to both dissipative and non-dissipative currents in the 
presence of weak magnetic field. Recently, some related 
studies have tackled this problem inspired by newly dis- 
covered Dirac electron materials, such as graphene and 
a-ET2l3, which are two-dimensional gapless Dirac elec- 
tron systems [Sl |9|, llOf. Although the anomalous behav- 
iors of X and the Hall coefficient Ru have been reported, 
their relationship has not yet been clarified. Here the 
relationship between x ^-^^id i?H in weak magnetic fields 
resulting from interband effects will be investigated, in 
detail, in the presence of the finite band gap 2 A. 

Several experiments on Rn of Bi have suggested that 



is quite anomalous, e.g., rapid sign change s If ll . 
141 ■ and 3D fractional quantum Hall effect [l5(. 



Rh 

spite of these fascinating behaviors, there has been no 
theoretical analysis especially in the weak field limit. For 
a two band system with a small band gap 2 A, such as 
Bi, it is puzzling how i?H behaves in the band gap region 
|/i| < A. Here, /x the chemical potential, and the origin 
of the energy is taken at the center of the band gap. In 
the free electron system, i?H is proportional to l/rioec. 
Tic being the density of electrons. Thus it should diverge 
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FIG. 1: Band structure of pure bismuth near L-point in the 
Brillouin zone. The chemical potential is /i ~ 27meV from 
the band-edge of the conduction band, and the band gap is 
2 A ~ 14meV01. 



for |/x| < A, since rio = 0. More precisely, if we start 
from the general expression, i?H = o'xy/HaxxO'yy, valid 
even in the multi-band cases, we obtain i?H ~* 0/0 for 
1^1 < A. (Here cr^j, is the conductivity tensor.) Namely, 
it is not trivial whether Rn converges or diverges and 
then the property of Rh in the band gap region is quite 
puzzling. 

In this Letter, we calculate axx, <^xy and i?H in the 
three-dimensional Dirac electron systems on the basis 
of gauge invariant Kubo formula in the Luttinger-Kohn 
representation [If 



17, igl. We find that the intcrband ef- 



fect induces unconventional contribution to dxy near the 
band-edge (i.e., |/i| '^ A), and i?H exhibits unexpected 
features; a sharp peak structure at around /i = ±A and 
a drastic sign change through n = Q. We discuss impli- 
cations of the present results to Bi and its alloys. 

The fact that low-energy properties of Bi can be de- 
scribed by Dirac electrons was first introduced by Co- 
hen and Blount [l9|, and then by Wolff in an elegant 
way[20J. For pure Bi, electrons locate near the L-point in 
the Brillouin zone and holes near the T-point. Previous 
studies concentrated on the L-point electrons, since their 
contribution is dominant due to much larger velocity of 
electrons than that of holes at the T-point. Based on 
this model, anomalous behaviors of Bi, such as a large 
g-factor[19J, magneto-optics [20j, large diamagnetism[J| 
have been clarified. Electrons around L-points (also holes 
around T-points) are essentially described by the tilted 
Dirac equations, i.e., Dirac equations with anisotropic 
velocitiesjsl, l20|. However in order to avoid complexities, 
which are not essential to present studies, the velocity of 
electrons, v, are assumed to be isotropic for simplicity. 
Then the Hamiltonian for electrons and holes is given by 



jr = A/3 -h iw V fcp/Ja^ 



(1) 



where (5 and ai are the 4x4 matrices that appear in the 
Dirac theory 
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with (Ti's being the Pauli spin matrices. This 
"relativistic" Hamiltonian already includes spin-orbit 
interactions [J, |20|, which is very strong in Bi contrary 
to graphene or a-ET2l3. 

It should be emphasized here that this Hamiltonian 
is not in the Bloch representation but in the Luttingcr- 
Kohn (LK) representation [l6|, in which the wave func- 
tion, -^nfc, is expressed in the form -0„fc = Unko{''')e^^'^ , 
Unko being the periodic part of the Bloch function at a 
wave number ko, e.g., the L-point for electrons or T-point 
for holes. The wave function under the magnetic field 
is correctly represented only in this representation 17 1 . 



Moreover, in the LK representation, we can obtain 
straightforwardly the gauge-invariant results for various 
correlation functions 171, ll8[ . Note that this represen- 
tation is exact and related to the ordinary Bloch wave 
function by a unitary transformation. 

The results of the conductivity, cFxx , smd the Hall con- 
ductivity, a^y, (for detail calculations, see Ref. [17|) are 
summarized as follows: 

2 l-oo poo 

ffx. - -^r / ^^/'(^ - ^') / ^^ [^1 (^' ^) - ^2 (e, X)] . 

TT-^W J-oo Jo 

(3) 



Ti(£,X) 



^2(^2^p2_ 1^2_^2) 



{(£ + ir)2 - X2 - A2} {(£ - ir)2 - X2 - A2} ' 

(4) 



, ^ X2{ £-fir2_lx2-A2} 

T2 e, X) = ^^ 2 i + c.c, 

^ ' 2{(£ + ir)2-X2-A2}2 



(5) 
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(6) 

(7) 



{(£ + ir)2-A2f/' 

-2r4 - r2A2 + (A2 - £2)2 + 2ir3£ - irM(A2 - £2) 



2r2£2 V£2 - r2 - A2 + 2ir£ 



(8) 



where f{e) is the Fermi distribution function. Here we 
have introduced a finite damping, F, for electrons as 
in Ref. [8|] to represent the effects of impurity scat- 
tering present in actual materials. (We assume F to 
be constant in order to make our argument as simple 
and transparent as possible, although Shon and Ando 
have indicated that F somewhat depends on energy and 
momcntum[2l|.) The chemical potential dependence of 
axx and axy at T = are shown in Fig. [2] (a) and (b) , 
respectively. The normalization factors are taken to be 
cTxxO = e^ JTV^v and Oxyo = e^v/12Tr^c. Away from the 
(.^) band-edge, |/i| » A, axx oc fi'^ and axy oc — ^, which are 
consistent with the Bloch band picture. In the band gap 
region, |^| < A, both conductivities have only small val- 
ues due to F. These non-zero values produce an anoma- 
lous i?H as shown below. 

The Hall coefficient defined by Rn — axy/a^^H is 
shown in Fig. [2] (d) where the normalization is taken 
to be i?Ho = '^xyo/'^xxo^^- Away from the band-edge. 



i?H oc ~ii^ as is expected from the properties of axx 
and axy This is consistent with the Bloch band picture 
since carrier density is proportional to \fjP\ in three di- 
mension. For \fj,\ < A, on the other hand, i?H exhibits 
unexpected features; a peak structure at ^ ~ ±A and 
a rapid sign-change through /i = 0. Note that there is 
no carriers at T = of |/i| < A in clean systems. It is 




FIG. 2: The chemical-potential dependences of (a) the con- 
ductivity a XX, (b) the Hall conductivity (Jxy, (c) the inter- 
band contribution to (j]^y''^ , and (d) the Hall coefficient _Rh at 
T = 0. The inset of (c) is that of the orbital susceptibility x. 



rather surprising that i?H docs not diverge but has small 
values (esp. -Rh = at /i = 0) for such small carrier 
density. The very particular feature, i?H = at ^ = 0, 
is similar to the case of gapless Dirac electrons in tvifo- 
dimensionj^, 0]. 

Now we study the interband effects of magnetic field 
on axy We extract the interband contribution, cr™''''', by 
subtracting the intraband contribution, ctJ,"*"^, from (Txy, 
namely, 



inter 

xy "x 



' xy 



(9) 



intraband approximation, i.e., the Bloch band picture (cf. 
[17|). which is given by 
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(10) 



where E±{X) = ±VA^ + A^. In spite of the appar- 
ently different expression between Eqs. ([6]) and (jlOp. they 
agree with each other except for the band-edge region; 
CT™*'''^ has only small value for |/i| < A, and increases as 



xy 
^intra 



-/z for |/i| > A. 



The obtained interband contribution cr^"*''' is shown in 
Fig. [2] (c) . The remarkable property is that crjf^*'"' takes 
the largest value at the band-edge, and becomes smaller 
away from the band-edge, contrary to ctJ,"*'''. Further- 
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Here cr^y''^ is the Hall conductivity calculated within the 



more, (J^y"" does not depend on i so much, while cr^^y^'^ 
does. (Note that the vertical axis of Fig. [2] (b) includes a 
factor r^.). This indicates that cr™*'''' has a different na- 
ture from cr^xy^^- ^^ seen in the inset of Fig. [2] (c), which 
shows orbital susceptibility x calculated in the gauge- 
invariant formula[l8J, cr^'L''"^ appears to be associated with 
diamagnetic currents. 

In the following, we discuss the physical origin of the 
interband contribution a'"*"'' in relation to the orbital 

xy 

susceptibility In the band gap region, electrons in a 
magnetic field circulate locally and make diamagnetic 
current. There are no electrons going through the crys- 
tal because of the band-gap. In the band-edge region, in 
contrast, the electrons start to be transferred from local 
diamagnetic orbitals to the neighboring orbitals due to 
some scattering. These processes will generate contribu- 
tions to the Hall conductivity: <7™y°^ ■ As is known from 
the behavior of x (Fig. E] (c)), the diamagnetic current 
has largest values for |^| < A, and decreases away from 
the band-edge |/i| » A. Correspondingly, |ct^'^*®''| actu- 
ally decreases away from the band-edge. Note that in the 
limit of A/r -^ 0, cr^"^'"' is given only by F^ term of Eq. 
([6|), where the factor f[e — y) generates the contributions 
of bands below //. 

Finally, let us discuss the implications of the present 
results to Bi. Some measurements of i?H in Bi have al- 
ready exhibited the peak structures and the rapid sign 
changes with respect to ^[1^, [14 1. They are, however, 
not consistent with each other. This would be due to the 
purity of samples or the inhomogeneity of either chem- 
ical or external pressure. For further investigations, we 
need some guidelines. It will be rather difficult to find 
the peak of Ra at // ~ ±A when the peak is very sharp. 
So we propose a simultaneous measurement of x together 
with that of i?H to observe the band-edge property in a 
more transparent way. A clear kink of y at /i ~ ±A has 
been actually observed experimentally [5J. Thus, the kink 
in X can be used to identify the band-edge, and then we 
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(c) ne > nh (electron doped) 



FIG. 3: Schematic band structure around L- and P-points of 
Bi and its alloys for (a) undoped (tic = nh), (b) hole doped 
{jic < nh), and (c) electron doped (ne > nh) case. The dashed 
lines indicate the level of the chemical potential. 



can expect the peak of Rh in the same region. 

When Bi alloys are used, it will be important to con- 
sider the contributions of both electrons and holes. The 
band structure around L- and T-points are illustrated in 
Fig. [3] for (a) undoped Bi (rie = rih), (b) hole doped 
Bi (rio < rih), and electron doped Bi (rie > nh). Here, 
X denotes some controlling parameter which changes the 
relative position of electron- and hole-band, for example, 
external pressure or alloy concentration. For undoped 
case (Fig. [3] (a)), a clear band-edge property — a peak 
in i?H and a kink in % — can be seen, since ric and rih 
vanish simultaneously at a; = xi. This band-edge prop- 
erty will be clearer when H is applied perpendicular to 
z (trigonal axis), since the contribution of electrons is 
dominant [J]. In this case, the sign change in Rn will 
not be observed, since n is always positive for electrons. 
For the hole doped case (Fig. [3](b)), on the other hand, 
/i of electrons can be negative for x > X2, so that the 
sign change in Ru is possible. However, the band-edge 
property in Ra will not be so clear due to the finite con- 
tribution from the holes at T-point, even though it will 
be clear in x- For the electron doped case (Fig. [3] (c)), 
neither the band-edge effect nor the sign change will be 
hardly seen. 

We have studied the Hall effects of the Dirac elec- 
tron systems with special reference to Bi on the basis of 
the Kubo formula in the Luttinger-Kohn representation 
keeping gauge invariance. It is clarified that the inter- 
band effect of a magnetic field generates the novel contri- 
bution to ax:y This contribution due to interband effect. 



cr™'*''', is remarkable near the band-edge. This property 
is quite different from the conventional Hall conductivity 
which originates from a single Bloch band. The Hall coef- 
ficient exhibits two unexpected properties; a sharp peak 
at // = ±A and a rapid sign-change through fi = 0. We 
can expect a clear structure of i?H even in the band gap 
region |^| < A. In order to detect the band-edge effect 
clearly, we have proposed simultaneous measurements of 
Ru together with x- We have predicted that the band- 
edge property and the rapid sign-change can be seen for 
the undoped Bi alloys (ue = ^h) and for the hole doped 
Bi alloys (ric < ^h), respectively. 
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